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We analyze the stability of the spontaneously dimerized phase of the frustrated Heisenberg anti- 
ferromagnet - the Ji — J2 model, in two dimensions. The lowest triplet excitation, corresponding 
to breaking of a singlet bond, is found to be stable in the region 0.38 < J2/J1 ^ 0.62. In addition 
we find a stable low-energy collective singlet mode, which reflects the spontaneous violation of the 
discrete symmetry. The spontaneous dimerization vanishes at the point of second-order quantum 
transition into the Neel ordered phase (( J2/ Ji)ci ~ 0.38). We argue that the disappearance of dimer 
order is related to the vanishing of the singlet energy gap at the transition point. 
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The nature of excitations in the quantum disordered 
phases of low-dimensional quantum antiferromagnets is a 
topic of fundamental importance for the physics of quan- 
tum magnetism [Q. Such phases can result from mobile 
holes in an antiferromagnetic background as in the t — J 
or Hubbard model at finite doping. Alternatively, com- 
petition of purely magnetic interactions can also lead to 
destruction of long-range order. A typical example of 
the second kind is the J\ — J2 model which exhibits a 
quantum disordered (spin-liquid) phase due to second- 
neighbor frustrating interactions. Even though it has 
been intensively studied during the last ten years, the 
J 1 — J 2 model still holds many secrets, especially concern- 
ing its quantum disordered phase. Exact diagonalization 
studies Q have shown that the excitation spectrum of 
the model is quite complex and that finite-size effects are 
large Spin- wave like expansions around the ordered 
phase (which occurs for small frustration) naturally can 
not give any information about the disordered phase at 
stronger frustration, and consequently non-perturbative 
methods are needed to analyze the latter regime. 

A new insight into the disordered regime was achieved 
by field-theory methods @J3], dimer series expansions 
j| and related effective Hamiltonian approaches HQ. 
The above works have shown that the ground state in 
the disordered regime is dominated by short-range sin- 
glet (dimer) formation in a given pattern (see Fig.l.). 
The stability of such a configuration implies that the lat- 
tice symmetry is spontaneously broken and the ground 
state is four-fold degenerate. Such a route towards quan- 
tum disorder is known rigorously to take place in one di- 
mension, where the Lieb-Schultz-Mattis (LSM) theorem 
guarantees that a gapped phase always breaks the trans- 
lational symmetry and is doubly degenerate, whereas 
gapless excitations correspond to a unique ground state 
||. A generalization of the LSM theorem to higher di- 
mensions has been suggested |J, but can not be rigor- 
ously proven. 

In this Letter we study further the stability of the spon- 
taneously dimerized phase in the J± — J 2 model and find 



that in the disordered phase not only the one-particle 
triplet excitations but also the two-particle singlet modes 
are stable. Both types of excitations are gapped and be- 
come gapless at the transition point to the Neel phase. 
The singlet collective mode is a low-energy excitation of 
the system and we argue that it is responsible for the 
decrease and ultimately the vanishing of the dimer or- 
der parameter as the transition is approached from the 
disordered side. Since the excitation spectrum is com- 
plex and close to the transition point the triplet is not 
the only low-energy excitation, it is quite possible, in our 
view, that the 0(3) non-linear sigma model in its usual 
formulation ]jj is not the correct effective field theory for 
the frustrated antifcrromagnet. 

The Hamiltonian of the Ji — J2 model reads: 

H = J 1 ^2s i .S J + J 3 ^Si.S j , (1) 

nn nnn 

where J\ is the nearest-neighbor, and J2 is the frustrating 
next-nearest-neighbor Heisenberg exchange on a square 
lattice (see Fig.l). Both couplings are antiferromagnetic, 
i.e. J\ t 2 > and the spins Si = 1/2. In order to study 
the excitations and their stability in the dimer phase we 
first derive an effective Hamiltonian in terms of bosonic 
operators creating triplets t\ a ,a = x,y,z from the sin- 
glets formed by each pair of spins, as shown in Fig.l. 
Similar effective theories have been derived in Refs. ||j7]] 
and we only present the result: H = H2 + H3 + H4 + Hjj, 

#2 = £ { Ma*k« + f 1 + } . ( 2 ) 

H 3= R ( kl ' k2 ) e «/3-r*k 1 a*k 2 /3*k37 + h ' C -> ( 3 ) 

1+2=3 

H 4 = T ( k i - k s) x ( 4 ) 

1+2=3+4 

(SasS/3-f ~ ^a/3 < V)*k lQ 4; 2 /3 i k37 i k4<5- 

We also introduce an infinite repulsion on each site 
{Hi/), in order to enforce the kinematic constraint on 
the Hilbert space t\ a t\g = 0. 
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The constraint is necessary in order to ensure that the 
bosonic Hamiltonian in terms of the triplet operators cor- 
responds uniquely to the original spin Hamiltonian (|l|) 
and no unphysical states appear in the final result. The 
introduction of Hjj into the effective theory and its sub- 
sequent diagrammatic treatment follow closely our pre- 
vious work Ref. Jl0| . 

The following definitions are used in Eqs.(||-||) 

Ak = J\ — ^-£fc x + {Ji - J2)£,k y - ■hS.k^ky, 

B k = - + ( Ji - J 2 )&, - M k J ky ■ (6) 

where = cos(fc), 7^ = sin(fc). The matrix elements in 
the quartic and cubic interaction terms are: 



T(k) = 



J1 + J2 



■h 



R(p,q) 



Jl 

~T 7px 



2 2 



(7) 



Throughout the paper we work in the Brillouin zone of 
the dimerized lattice. 

The spectrum is determined by the normal Green's 
function (l(J: G N {k,oj) = [u + A k (-w)][{w + 
A k (-uj)}{uj - A k (w)} + Bl]- 1 where i k (w) = A k + 
Xjv(k, oj) and i?k = £>k + Syi(k, w). The normal self- 
energy Ejv is given by the diagrams of Fig. 2(a) where 
the effective scattering vertex T has to be found by solv- 
ing the Bethe-Salpeter equation shown in the same fig- 
ure. The first two diagrams in Fig. 2(a) (proportional to 
r and T) can be calculated as explained in Ref. p0[ , 
while the three-particle contribution (third diagram in 
Fig. 2 (a)) was discussed in more detail in Ref. jl]]] . We 
adopt the Brueckner approximation which is based on 
the smallness of the triplet density n t = ^2 a {t\ a ti a ) -C 1 
(Ref. jnj). However for the system under considera- 
tion n t ~ 0.3 is not very small and additional dia- 
grams have to be summed up to ensure the reliability 
of the results. This mainly concerns the anomalous self- 
energy Y,a shown in Fig. 2(b). One improvement was dis- 
cussed in Ref. 12 and involves summation of rainbow 
diagrams containing U . This is equivalent to introduc- 
tion of the term AJ~j. (t] a t| Q + h.c.) into the Hamilto- 
nian and choosing the Lagrange multiplier A from the 
condition {t\ a t\ a ) = 0. The second improvement is the 
summation of rescattering diagrams which take into ac- 
count the attraction in the singlet channel (see Fig. 2(c) 
and the discussion of the singlet bound state below). The 
two effects put together give 
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(5) The pole of Gff(k.,u), which has to be found self- 
consistently, determines the renormalized one-particle 
spectrum w(k). Here we only present the results, ob- 
tained by a numerical iterative procedure. The dot- 
dashed line in Fig. 3. is the spectrum for the particu- 
lar value of J2/J1 = 0.4. For comparison we have also 
performed a dimer series expansion up to order 8, and 
have resummed the dimer series by using Pade approxi- 
mants |l3| . The result for the triplet excitation spectrum 
is plotted as a solid line in Fig. 3. The spectra obtained 
by the two methods are in excellent agreement. One can 
also see that in some regions of k space w(k) lies in the 
shaded region in Fig. 3, which represents the two-particle 
scattering continuum. The excitations can decay in that 
region which is reflected in the poorer convergence of the 
dimer series (larger error bars on the solid curve). 

To map out the phase diagram of the model we have 
studied the instability of the spectrum as J%/J\ varies. 
At the critical value (J2/Ji) c i = 0.38 the gap at k = 
(0, 7r) (which corresponds to the Neel ordering wave- 
vector (71", 7t) of the original (non-dimerized) lattice) van- 
ishes, signaling a transition to the Neel ordered phase. 
For larger frustration the gap grows almost linearly with 
J2/J1, as shown in Fig. 4(a). At (Ji/Ji)c% = 0.62 an in- 
stability towards a collinear phase takes place and the 
gap at (0, 0) ((tt, 0) of the original lattice) goes to zero. 
The transition at (J2/ Ji) C 2 appears to be of first order 
and is quite similar to the transition which occurs in the 
Heisenberg ladder with frustration The number of 
triplets increases sharply near this transition point and 
the excitation spectrum is a complex mixture of one- and 
many-triplet bound states |L^ |. In what follows we will 
concentrate on the vicinity of the second-order transition 
into the Neel phase. 

Let us now investigate the spectrum of collective ex- 
citations - in particular we have studied the two-particle 
bound state with S=0. This state has the form: |\&q) = 

E q *(q:Q)^,Q/2+ q t i,Q/2- q l°}' where Q and q are > re- 
spectively, the total and relative momenta of the two 
quasiparticles. The wave function ^ satisfies the inte- 
gral equation (£' 5 (Q) is the bound state energy) JTof : 

[£ S (Q) - w Q/2+q - uj Q /2-q\ *(q, Q) - 

{-2[T(p - q) + T(p + q)] + U}#(p, Q). (9) 

p 

Since we have to take U — > 00, the following con- 
dition must be imposed via a Lagrange multiplier: 
^ q \E'(q, Q) = 0. The numerical solution of Eq.(||) is 
presented in Fig. 3 with a long dashed line. Notice that 
the singlet bound state has a very low energy and exists 
below the two-particle continuum for all wave-vectors. 
The singlet is gapped everywhere in the disordered phase, 
however its energy at k = (0, 0) approaches zero near the 

&.&,,+ A. (8) 

transition point to the Neel phase as shown in Fig. 4(a). 
We emphasize that the appearance of a low energy sin- 
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glet in the spectrum is quite unusual and represents a 
characteristic feature of frustrated spin systems. 

Next we calculate the dimer order parameters in the 
two spatial directions, defined as (see Fig.l for notations): 
D x = (S2.S3) - (Sx-Sa), D y = (Sx.Sg) - (Si.S a ). Dimer 
series expansions to order 9 have been carried out for 
these quantities and the results are presented in Fig. 4(b). 
Both order parameters appear to approach zero at the 
transition point to the Neel phase, even though the error 
bars are also quite large. Our result is very different form 
earlier works || where no substantial decrease of the or- 
der parameter was found. We attribute the difference to 
the longer series we have generated. We have also calcu- 
lated the dimer order parameter by using the diagram- 
matic approach. The results shown in Fig. 4(b) suggest 
that the dimer order is overestimated in this way since 
it shows no appreciable decrease as frustration decreases. 
However we envisage the following physical mechanism 
which would explain the discrepancy between the two 
calculations. As discussed earlier, when the transition 
point (J2I Ji)ci = 0.38 is approached, the energy of the 
singlet state decreases and ultimately goes to zero. Since 
the ground state is also a singlet, a strong mixing between 
the bound state singlet and the ground state must occur 
near the transition, which effectively leads to the increase 
of quantum fluctuations. Technically this effect can be 
taken into account by exact calculation of the diagram 
shown in Fig. 2(b). This diagram represents a contribu- 
tion to the anomalous Green's function and thus effec- 
tively to the strength of quantum fluctuations. The ex- 
act two-particle scattering amplitude M which appears in 
Fig. 2(b) satisfies the equation shown in Fig. 2(c). Notice 
that the bound state equation (|^) is the pole of M, and 
it is quite easy to prove that as the bound state energy 
decreases, the contribution of the diagram in Fig. 2(b) in- 
creases substantially. The numerical implementation of 
this procedure is however quite difficult. The renormal- 
ization of given by Eq.(|8j) is a step in this direction, 
but does not fully take into account the effect. Never- 
theless it is clear that the vanishing of the singlet gap 
at the transition is intimately related to the increased 
quantum fluctuations and therefore the vanishing of the 
dimer order parameter. 

Finally we present exact diagonalization data on a 
small cluster which provide further evidence that the 
ground state is spontaneously dimerized. Since there are 
four ways the system can choose the dimer pattern in 
Fig.f, the ground state is expected to be four- fold de- 
generate. Indeed, three singlet states are clearly seen 
above the ground state (See Tabled) even though they 
have non-zero excitation energy on a finite cluster (an 
exponentially small splitting is expected between them 
for large cluster size). One also expects that upon in- 
troduction of explicit dimerization 5 into the system (i.e. 
coupling (1+5) on the bond [12], [34], etc. (see Fig.l.) 
and (1-5) on the bonds [23], [15], etc.), the degeneracy 



of the ground state will be lifted and therefore the three 
zero energy singlets, related to the spontaneous dimeriza- 
tion, should disappear. On a finite cluster the energies 
of these singlets should grow as 5 x (system size) . Indeed 
one can see from Tabled that the singlet energies increase 
with 5, which proves their symmetry breaking origin. 

In summary, we have explicitly calculated, for the first 
time, the lowest triplet and collective singlet excitations 
of the two-dimensional J\ — J2 model. Both branches 
of the spectrum were found to be stable and gapped 
in the quantum disordered phase and become gapless at 
the boundary with the Neel ordered phase. We empha- 
size that a low-lying singlet around k = (0, 0) is present 
only in a spontaneously dimerized system. In contrast, 
for systems with explicit, exchange driven dimerization, 
the binding energy for such a singlet it strictly zero. 
Thus the low-energy singlet around the Brillouin zone 
center should be, in principle, observable in quasi two- 
dimensional systems exhibiting a quantum critical point, 
separating a dimerized (spin-phonon interaction driven) 
and a Neel phase [O]. In addition, we have found that 
the dimer order parameter tends to zero at the transi- 
tion point and have argued that the disappearance of 
spontaneous dimer order is tied to the vanishing of the 
singlet gap. Let us mention that the structure of the 
spectrum, found in the present work, supports the pic- 
ture of Read and Sachdev [Q , based on their analysis of 
the Sp(N) field theory. In particular, both in our and in 
the field-theory approach, a singlet mode associated with 
the spontaneous dimer order is present. 

The example of the J\ — Ji model provides further sup- 
port for the idea, put forward some time ago || , that the 
Lieb-Schultz-Mattis theorem can be extended to higher 
dimensions, and the gapped states of quantum systems 
necessarily break the discrete symmetries of the lattice. 
In this connection it would be very interesting to inves- 
tigate if similar conclusions can be reached for systems 
with finite doping, e.g. the t — J model away from half 
filling. The recent discovery of stripes in the high-T c ma- 
terials Jl^], as well as the analysis of the Sp(2N) t — J 
model for large TV strongly suggest that the disor- 
dered ground states at finite doping may be unstable to- 
wards configurations which break the lattice symmetries. 
Further analysis of these possibilities is highly desirable. 
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TABLE I. Exact diagonalization results on N=4 x 4 cluster 
for the energies of the lowest five excited states at k = (0, 0). 
Frustration is fixed at J2/J1 = 0.5 and S is the explicit dimer- 
ization. "s" and "t" denote singlet and triplet levels. 



5 = 



0.319 (s) 
0.830 (s) 
0.971 (s) 
1.128 (t) 
1.135 (s) 



0.1 



0.655 (s) 
0.864 (s) 
1.001 (t) 
1.144 (s) 
1.425 (t) 



0.2 



0.908 (t) 
0.924 (s) 
1.327 (s) 
1.506 (s) 
1.684 (t) 



0.3 



0.864 (t) 
1.130 (s) 
1.814 (s) 
2.001 (s) 
2.028 (t) 



FIG. 1. Schematic representation of the system. The ovals 
represent two spins coupled into a singlet. 



FIG. 2. (a) Lowest order diagrams for the normal 
self-energy Sjv. (b) Contribution to the anomalous 
self-energy Ea- (c) The Bethe-Salpeter equation for the 
two-particle scattering amplitude M. 



FIG. 3. Excitation spectrum for J2 = 0.4Ji. The 
dot-dashed and solid lines are the diagrammatic and dimer 
series results for the triplet dispersion, respectively. The long 
dashed line is the two-particle singlet excitation. The shaded 
area represents the two-particle scattering continuum. 



FIG. 4. (a) Triplet and singlet energy gaps in the disor- 
dered phase. Solid lines represent the dimer series results for 
the triplet gaps. The dotted and dashed lines are the diagram- 
matic calculations of the singlet and triplet gaps, respectively, 
(b) The dimer order parameters D x and D y (see text for defi- 
nitions). Dashed lines are the diagrammatic results and solid 
lines are the dimer series results. 
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Fig.l. 
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Fig.2. 




(k x , k y ) 



Fig. 3. 



